Self-Polarization and Cooling of Spins in Quantum Dots 
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Spontaneous nuclear polarization is predicted in double quantum dots in the spin-blocked electron 
transport regime. The polarization results from an instability of the zero-polarization state when 
singlet and triplet electron energy levels are brought into resonance by the effective hyperfine field 
of the nuclei on the electrons. The nuclear spins, once polarized, serve as a cold bath for cooling 
electrons below the lattice (phonon) temperature. We estimate the relevant time scales and discuss 
the conditions necessary to observe these phenomena. 



Recent advances in semiconductor quantum dot tech- 
nology have given experimentalists the ability to con- 
trol the behavior of individual electrons and investigate 
their coupling to nuclear spins. Interesting phenomena 
such as switching, hysteresis, and long period oscillatory 
behavior of electric current were observed by Ono and 
TaruchaQ in the so-called spin-blockade regime in GaAs 
vertical double quantum dots. Koppens et al. have also 
observed bistability and switching of electron current in 
lateral quantum dot system under analogous conditions 
0. In both cases, strong evidence was presented link- 
ing the observed phenomena to collective behavior of the 
nuclear spins in the lattice. 

Collective behavior of nuclear spins in solids at Kelvin 
temperatures [li, suggestive of spontaneous ordering, 
is a remarkable phenomenon since the natural dipole- 
dipole interaction between nuclear spins is so weak that 
microkelvin temperatures are required to achieve order- 
ing under equilibrium conditions. The self-polarization 
of nuclei observed in Refs.P,|3] was made possible by the 
non-equilibrium conditions of these experiments. 

Much of recent theoretical work on spins in quantum 
dots has focused on dephasing fsl ^ and relaxation 
of electron spins due to their interaction with the lattice 
nuclei. Additionally, Jouravlev and Nazarov^ proposed 
a theory of spin-blockaded electron transport governed 
by the effective field of nuclear spins. In contrast to the 
idea of random fluctuations of the disordered nuclear spin 
system, in this work we are interested in the situation 
where spin-blockade leads to spontaneous ordering of a 
macroscopic fraction of the nuclear spin system. 

Although qualitative scenarios explaining the above- 
mentioned instabilities have been offered^ Q, a com- 
plete picture remains elusive. In a different context, 
however, Dyakonov and Perel studied instabilities involv- 
ing nuclear spins in the optical pumping regime in bulk 
semiconductors • Between this theoretical work and the 
unexplained observations of current-driven instabilities, 
there is a gap to be filled in. 

In this paper we discuss a model of electron transport 
through a double quantum dot system which exhibits an 
instability toward self-polarization of the nuclei. In equi- 
Ibrium at temperatures above a few millikelvin, nuclear 
spins are randomly oriented, having a root-mean-squared 



polarization of size ©(v N), where is the total number 
of spins. To achieve an extensive polarization that scales 
with A^, the system must be driven out of equilibrium. 

Polarization of nuclear spins occurs as an unpolarized 
electron current is passed through the double dot device. 
The hyperfine interaction between electrons and nuclei 
couples the dynamics of electron and nuclear spins. In a 
certain window of level detuning, feedback of the effec- 
tive field of nuclear spins on electron levels near resonance 
becomes strong enough to pull the levels further into res- 
onance by polarizing the nuclear spins. Transport drives 
the instability, taking the nuclear spins out of equilib- 
rium and allowing an extensive polarization to develop. 
We discuss this effect and show that it can give rise to 
collective behavior of the otherwise independent nuclei. 

Once one has successfully polarized the nuclear spins 
within a quantum dot, an interesting thermodynamic sit- 
uation is obtained. Consider the equilibrium of a sin- 
gle electron spin isolated in such a quantum dot in the 
presence of an applied Zeeman field. While phonons 
are in equilibrium at the lattice temperature, the (non- 
equilibrium) nuclear subsystem has an effective tem- 
perature that is much lower. Assuming a static, non- 
equilibrium nuclear spin background, to what effective 
temperature does the electron spin relax? 

The non-equilibrium nuclear spin distribution can af- 
fect the electron spin dramatically if electron spin re- 
laxation is predominantly due to hyperfine spin flips. 
We compare the contributions to relaxation due to spin- 
orbital and hyperfine couplings, and conclude that the 
triplet states of a double dot system with two electrons 
can be cooled to an effective temperature Tegf below that 
of the phonon bath [see Eq. IjlUll ]. 

In the paragraphs that follow, we outline the model put 
forth by Ono and Tarucha in [ij to explain the observed 
behavior of their vertical double quantum dot devices 
in the spin blockade regime. This model, accepted here 
for concreteness and to set the stage for our discussion, 
illustrates more general ideas that apply to a variety of 
systems with similar electron energy spectra. 

The double dot is weakly coupled in series to two un- 
polarized leads, with transport occuring as a series of dis- 



crete hopping events as described in Figure 1 (a) Current 
is suppressed when one of the three (1, l)t states are oc- 
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Fermi's Golden Rule: 
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(b) Relevant Transitions 



FIG. 1: (a) Schematic of the double quantum dot system. 
Initially the system is in state (0, 1) with one electron on 
the right dot. Current flows as depicted by the arrows: an 
electron tunnels from the source into the left dot to form the 
state (1, l)s/t with singlet or triplet spin, respectively. The 
electron then must hop to the right dot to form (0, 2)s before 
tunneling out and returning the system to the (0, 1) state. If 
(1, l)t is formed when the left dot is filled, the Pauli Principle 
prohibits the second electron from tunneling into the right 
dot. (b) The relevant transitions included in our calculation. 
Bare spin flip rates W± are calculated by Eq.Q. Electrons 
can also escape from the dot without nuclear spin flip with 
rate PF'", giving the net spin flip rate of the form Q. 

cupied in the first step of the cycle. Once in the triplet 
state, residual current is due to slow indirect tunneling 
through virtual excited states, exchange with the leads, 
and spin-flips due to spin-orbit coupling and / or hyperfine 
flip-flop scattering. However, because spin-orbit effects 
are suppressed due to confinement in structures of this 
type Q> Q > we consider only the effects of indirect 
tunneling and hyperfine scattering. 

For simplicity, we assume that tunneling rates from the 
weakly coupled source lead into the four (1, l)^/^ states 
are all equal. Coupling of the (0, 2)s state to the drain 
lead is stronger and gives this level a significant decay 
width. The finite width of this level plays a key role 
in the feedback mechanism that drives the polarization 
instability by favoring transitions from one of the split- 
triplet levels as it approaches resonance with (0,2)^. 

In general, the orbital eigenstates with singlet spin con- 
figuration are a superposition of the states (1, l)s and 
(0, 2)s. Here we assume, without loss of generality, that 
one eigenstate retains predominantly the (1,1) charge 
distribution, while the other retains predominantly the 
(0,2) charge distribution. It is these "(1, l)s-like" and 
"(0, 2)s-like" states to which we refer below. 

The transitions relevant to our model are depicted in 
Figure 1(b) We assume that nuclear dynamics are in- 



coherent, and describe them by the populations N± of 
the up and down nuclear spin states [1^. We also ne- 
glect spatial variations in the nuclear spin population and 
transitions that do not change the net spin, which are 
inessential for our analysis. The energy-dependent hy- 
perfine spin fiip transition rates are calculated using 
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where /(e) is the density of states for the singlet final 
state, and et± , are the energies of the triplet and sin- 
glet, respectively. We assume a Lorentzian lineshape 

/(e) « 
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to allow explicit calculation. 

Our use of Fermi's Golden Rule is valid when nu- 
clear spins dephase on the time scale of a few electron 
tunneling events. This situation is to be contrasted with 
the picture involving coherent hyperfine coupling dis- 
cussed by Taylor et al. 16] which ignores nuclear spin 
dephasing and obtains typical polarization 0{^/N). 

When electrons are injected from an unpolarized 
source and every electron must exchange its spin with 
a nucleus to escape, no spin can be pumped into the nu- 
clear spin system irrespective of the ratio of the rates 
for flipping nuclei up or down, /W^ . If, on the other 
hand, electrons have an alternative way to escape, it need 
not be the case that the same number of nuclei must flip 
their spins in each direction. 

For simplicity, we assume a single energy-independent 
indirect tunneling rate W™ for all three triplet states 
that relieves spin-blockade without interaction with the 
nuclei. Due to the competition between these processes, 
the net nuclear spin flip rates are given by 
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where Iq is the total current through the system[l^. If 
the energy dependent rates are not equal, it is pos- 
sible for electron transport to be dominated by spin-flip 
processes for electrons of one spin type, and by indirect 
tunneling processes for electrons of the other type. 
The nuclear system is in a steady state when the rates 
and are equal. Assuming no dependence of the 
matrix elements containing the orbital wavefunctions on 
the electron spin-projection, we have (cf. Ref.0): 



/(£+)7V_ =/(e_)iV+ 



(3) 



When f{e+) ^ /(£-), there can be a nonzero nuclear 
spin polarization in the steady state even when electron 
Zeeman energy is negligible. Let x = (iV+ — N-)/N be 
the fractional nuclear polarization. Due to the hyperfine 
coupling-induced Overhauser shift, the triplet levels are 
split: £± = Eo ± aNx, where Eq is the detuning from the 
triplet-singlet resonance. For GaAs, where « —0.44, 
a < 0, this leads to the situation depicted in FigI21 

The condition (PJ yields a third order equation for the 
equilibrium polarization x* with solutions 
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FIG. 2: Electron energy levels involved in the spin blockade. 
The m = triplet level is detuned from the singlet by energy 
eo- When the nuclear spin system has fractional polarization 
x, the effective (Overhauser) field causes a splitting of aNx 
between the triplet levels. Resonance of one of the triplet 
levels with the broadened (0,2)s state favors spin flips for 
this state. 

Which of these solutions corresponds to a stable equi- 
librium? The solution a;* = always exists, for at 
-Bext = the triplet levels are degenerate when a; = 0. 
However, the solutions with finite nuclear polarization 
exist only when the discriminant is positive, i.e. 



(5) 



For a detuning of order 7, relation (O requires that 
the electron hyperfine shift at full nuclear polarization is 
larger than the singlet level width. When this condition 
is met, the zero-polarization solution x* = is unstable. 

Beginning at large detuning where is not satisfied, 
the system stably remains at the equilibrium x* — ^. 
This state loses stability once eo is reduced past the 
threshold, and the system tends to one of the solutions 
x^. At optimal detuning when eg = —oiN, the sponta- 
neous polarization reaches 



N+ - N-U^^ - ±^N^ - 72/^2. (6) 

This should be contrasted with an 0{'\fN) polarization 
due to the initially randomized nuclear spins. In typical 
GaAs quantum dots, fully polarized nuclei can create an 
effective field of strength of up to a few Tesla, which 
translates into a triplet splitting of et+ — £t_ ~ O.lmeV. 

For a total current Iq ~ IpA, there will be one spin- 
flip for every few tunneling events leading to a net spin- 
flip rate of 0{h/e) « 10^ - lO^s^i. In a sample with 
N — 10^ nuclei, this means that a particular nucleus is 
pumped with a rate 0{Io/eN) « 10 - lO^s-^ Nuclear 
relaxation is heavily suppressed at low temperatures, and 
can be assumed to take place on a much longer time- 
scale than this. However, the nuclear dipole-dipole flip 
rate Wdd ~ 10^s~^ (corresponding to a dipole field of 
1Q~^T) is non-negligible on the predicted time scale of 
pumping and must not be ignored. 

The dipole-dipole interaction does not conserve total 
spin due to the presence of terms such as S^S^ , S^S^ : 
etc. (e.g., see IHj). In a weak applied field, the Zee- 
man splitting of nuclear spin states with different total 
z-projection quenches the action of such terms due to en- 
ergy nonconservation. The strength of the applied field 



required for that need only give rise to Zeeman splitting 
greater than the dipole-dipole level width. Typically, a 
field Bcxt ~ 0.5mT is sufflcent for this purpose. 

The terms remaining after application of the small 
quenching field allow nuclear spin to diffuse throughout 
space (seeRef.[l3l for discussion of spin diffusion) . Using 
experimentally determined diffusion parameters |lll Il2l | , 
we estimate a characteristic time of 10 seconds for diffu- 
sion out of the dot. This rate is several orders of mag- 
nitude smaller than the rate of pumping and does not 
inhibit the pumping of nuclear spins. 

We now discuss the cooling of electron spins by the po- 
larized nuclear spin bath. Due to the small nuclear Zee- 
man energy compared to that of an electron, hyperfine- 
induced electron spin relaxation must involve additional 
degrees of freedom such as phonons to carry away the 
excess energy. Electron spin relaxation in a random nu- 
clear field in a single quantum dot has been studied by 
Erlingsson, Nazarov, and Fal'ko Spin-orbit relaxation 
rates were studied by Khaetskii and Nazarov ^ and by 
Golovach, Khaetskii, and Loss 0. These estimates along 
with recent experimental results |l3l | indicate that on 
time scales of interest the rate of spin relaxation for a 
single electron is dominated by the spin orbit mechanism 
with phonon emission. 

In contrast, Johnson et al. found that relaxation of a 
two-electron triplet state in a double quantum dot was 
dominated by hyperfine interactions at magnetic fields 
less than approximately 100 mT The reason for this 
reversal of roles is that hyperfine transitions are enhanced 
by the near resonance of singlet and triplet levels in the 
double dot system, whereas the spin orbit coupling is 
suppressed by a small tunneling matrix element. Thus in 
this low-field regime, the presence of a polarized nuclear 
spin bath may cause the electron spin degrees of freedom 
to adopt an equilibrium distribution different from the 
thermal distribution at the lattice temperature |2d| . 

We estimate the cooling effect for the system consid- 
ered above, with couplings to the leads tuned to zero and 
a gating bias favoring the unbroadened (0, 2)s state. As 
before, we treat transitions as incoherent using Fermi's 
Golden Rule, this time via second order perturbation 
theory in the hyperfine and electron-phonon couplings. 
At this level the transitions with Am = ±2 are not al- 
lowed. Rates for all other "forward" and "reverse" pro- 
cesses along with the normalization condition — 1 
yield a system of five linear equations in the five unknown 
equilbrium populations Pk, k = 0, 4 stands for (0, 2)s, 
(l,l)t+, (l,l)to,(l,l)t_, or (1,1)5, respectively. 

The rate equations governing the relaxation dynamics 
of this system are given by 



Pk = 



E 



(7) 



where Tkj is the transition rate from state j to state 
k. For example, the relaxation rate Fgi from (1, l)t , to 
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(0, 2)s via virtual transition through (1, l)s is 
Toi = rph[A£;oi] (M/AEiif 7V_, 



(8) 



derived as in Ref.j^ for the case of a single dot. Here 
AEij — Ei — Ej and rph[A£'] is the rate of phonon ab- 
sorption/emission for positive/negative energy AE. 

The phonon rates rph[A_E] in equation obey 
the detailed balance relation rp\i[AEjk]/Tp\i[AEkj] = 
g-Pph^Ejk ^ with Tph = the phonon bath tempera- 
ture. The ratios of forward and reverse rates for tran- 
sitions that include a change of m, however, include an 



additional factor of N^/N^ 



-/3„-E„ 



or its reciprocal. 



This relation defines the effective temperature T„ = /3„~ 
of the nuclear subsystem, with En = fJ.nB the Zeeman en- 
ergy of a single nucleus. Here /i„ is the nuclear magnetic 
moment and B is the applied magnetic field strength. 

In the steady state with Pk — for all k, the popula- 
tions obtained from equation (O obey 



P2/P0 = e-z^ph^-E^o 
P4/P0 = e-/3phA-B4o^ P3/P2 = e 



These relations describe a non-thermal distribution such 
that the energy dependence cannot be accounted for 
by a single effective temperature Tcff — PcS^- 

Non-Boltzmann aspects aside, this process can still be 
considered as a form of cooling as a finite (positive) nu- 
clear spin temperature leads to a transfer of population 
from high energy to low energy states. The energy lev- 
els within the triplet subspace are evenly spaced with 
AE21 — Ai?32 — fJ-eiB^ where fiei is the electron mag- 
netic moment. This leads to a quasi-thermal distribution 
within this subspace: e"'^""^-^"! ee P2/P1 = P3/P2 = 
g-0phAB2i g-i3„E„^ with effective temperature 
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(10) 



with fi„/lJ-ei ~ lO^'^. For a lattice temperature of 100 
mK and total field of 100 mT, Toff will differ significantly 
from T for nuclear polarizations greater than about 50%. 

To summarize, the key ingredient that makes it pos- 
sible to generate a finite nuclear polarization using cur- 
rent from unpolarized leads is the competition between 
energy dependent spin-fiip and energy independent spin- 
conserving tunneling channels. When spin-flips are fa- 
vored in one spin state but spin-conserving tunneling fa- 
vored for the other, polarization will be pumped into the 
nuclear spin system. Once a non-equilibrium nuclear spin 
state is established, it can be used as a cold bath to cool 
electron spins below the lattice temperature. While rela- 
tively fast spin-orbit processes in GaAs prohibit cooling 
of a single electron spin, two-electron states of a double 
quantum dot are amenable to such cooling. 

Electron transition dependence on Overhauser shift as 
a cause of bi-stability in spin-blockaded transport was 



discussed by Inarrea, Platero and MacDonald in Ref . j 
which appeared after this work was completed. 
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Although our primary example is GaAs which is com- 
prised of spin-3/2 nuclei, for clarity of presentation we 
proceed with the calculation assuming spin-1/2 nuclei. 
This allows us to introduce the populations of nuclei in 
the up and down spin states, N+ and N- . 
In principle the current can be found self-consistently 
from the rates of the individual processes. The energy 
and nuclear-spin state dependence of the current leads 
to a quantitative change in the rate of spin pumping, but 
does not affect our general conclusions. 
Although the Overhauser field of a nearly polarized nu- 
clear system is of order 1 T, an opposing external field 
can be applied to make the total effective field seen by 
the electrons fall into the low-field regime. 



